MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  OF  STANDARDS  A 


AD  No. 

DDC  FILE  COPY  ADA056488 


<g 


JIADIAL  AND  TRANSVERSE  RESPONSE  OF 
\ GUN  TUBES  TO  ^RAVELING’  BALLISTIC 
, PRESSURE  ; 


■cm 


/OJ 

— ■'7 


•78 


THOMAS  E.  SIMKINS,  RHfe  RAYMOND  D.  SCANLON,  MR. 

' GEORGE  A./PFLEGL/  1 

'"B^ETlVEXPO^STAKjRATORY 
WATERVLIET  ARSENAL,  WATERVLIET,  NY 


Introduction 


C - 


JVrJ 

DDC 

ElSEDIlCEfnl 

JUL  12  1978 


1SMTTE 


The  design  and  analysis  of  gun  tubes  presently  proceeds  on 
the  basis  of  the  static  equations  of  continuum  mechanics.  Stresses 
within  the  tube  wall,  for  example,  are  calculated  by  the  classical 
Lame  formula  - generally  viewed  as  a conservative  design  criterion. 
The  advent  of  high  speed  computation  and  mathematical  modeling,  how- 
ever, has  revealed  that  dynamic  bore  expansion  during  the  interior 
ballistic  cycle  will  create  significantly  higher  wall  stresses  than 
those  on  which  the  tube  design  is  based  The  analytical  and 

experimental  results  leading  to  this  conclusion  comprise  the  material 
in  the  first  section  of  this  work.^ 

The  second  portion  of  the  work  which  follows  shows  how 
transient  bending  vibrations  may  arise  during  firing  due  to  tube  cur- 
vature. Tube  motions  predicted  at  the  muz&le  are  of  sufficient  magni- 
tude to  explain  a portion  of  the  error  realized  at  the  target.  Three 
sources  of  curvature-induced  vibration  are  derived  (3)  and  two  are 
treated  in  detail. 


An  appreciable  effort  has  been  made  in  the  interest  of  real- 
ism. Highly  detailed  tube  geometries  and  interior  ballistic  curves  of 
pressure  and  projectile  travel  for  specific  weapons  of  interest  have 
been  included  in  the  analysis.  Wherever  possible,  use  has  been  made 
of  the  large,  widely  accepted,  NASTRAN  (NAsa  STRuctural  ANalysis) 
finite  element  computer  code.  Though  NASTRAN  is  quite  versatile,  it 
is  not  particularly  well  suited  for  handling  curvature -induced  load 
functions  which  require  special  programming. 
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1.  RADIAL  VIBRATIONS  IN  RESPONSE  TO  TRAVELING  BALLISTIC  PRESSURE 

The  175  mm,  M113  gun  tube  was  chosen  for  analysis.  The  tube 
is  assumed  to  be  completely  free  of  supports,  closed  at  the  breech  and 
subjected  to  a time-variant  internal  ballistic  pressure  which  traverses 
the  bore  at  projectile  velocity.  Of  primary  interest  is  the  radial 
displacement  of  the  bore  surface  of  the  tube  in  time.  To  this  end  a 
250  degree-of- freedom  model  of  the  tube  is  created  by  dividing  the  tube 
into  62  trapezoidal  ring  elements  as  shown  in  Figure  1.  These  finite 
elements  allow  the  true  geometry  of  the  tube  to  be  modeled  very  closely. 

Program  output  in  terms  of  radial  or  axial  displacements, 
velocities,  or  accelerations  may  be  requested  for  any  bore  station 
along  the  tube.  As  representative  output,  the  radial  displacement  vs. 
time  at  three  locations  along  the  tube  have  been  selected  for  presen- 
tation. In  all  there  are  125  points  in  the  model  at  which  program  out- 
put is  computed. 

Geometrical  information  for  the  NASTRAN  program  was  obtained 
from  tube  manufacturing  drawings.  Tabulations  of  the  pressure-travel 
curves  of  Figure  2 provided  the  input  necessary  for  load  specification 
(4).  Traveling  loads  cannot  be  prescribed  directly  for  NASTRAN  input 
and  require  special  programming  considerations. 


Figure  1 - Section  of  Trape-  Figure  2 - Ballistic  Curves  - 175  mm 
zoidal  Ring  Element.  M113,(Zone  3). 
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section  of  the  tube. 


RESULTS 

Figure  3 displays  the  NASTRAN-pre- 
dicted  radial  displacement  of  the 
bore  surface  at  a location  in  the 
breech  of  the  gun  tube.  Superposed 
for  comparison  is  the  displacement 
as  calculated  from  the  Lame  design 
formula  and  the  applied  pressure. 
Owing  to  its  location  to  the  rear 
of  the  projectile  at  all  times, 
this  station  receives  pressure 
smoothly  and  continuously  from  time 
zero  (ignition).  Consequently,  the 
‘ loading  rate  at  this  point  depends 
only  on  the  rate  at  which  the  burn- 
ing propellant  produces  pressure 
and  not  primarily  upon  the  velocity 
at  which  the  load  follows  the  pro- 
jectile. 


It  may  also  be  observed  from  figure  3 that  there  is  a stead- 
ily increasing  amount  of  natural  frequency  excitation.  These  vibra- 
tions have  a visible  tendency  to  form  'beats'  in  the  classical  sense 
and  infer  the  coexistence  of  several  modes  of  vibration  having  nearly 
equal  frequencies.  Their  importance  lies  in  that  the  amplitudes  of 
the  beat  envelopes  can  attain  quite  large  values,  creating  an  oppor- 
tunity for  overstress  and  shortened  fatigue  life. 

Finally,  it  is  noted  that  the  Lame-predicted  displacement 
response  is  conservative  in  that  its  maximum  value  exceeds  that  pre- 
dicted by  the  NASTRAN  analysis.  This  is  to  be  expected  since  the  Lame 
solution  is  formed  on  the  assumption  that  the  tube  is  uniformly  pres- 
surized throughout  its  entire  length.  In  reality  - and  in  the  NASTRAN 
model  - the  pressure  only  exists  in  the  region  to  the  rear  of  the  pro- 
jectile. The  unloaded  region  in  front  of  the  projectile  tends  to 
assume  its  unloaded  configuration  exerting  a contractile  effort  on  the 
pressurized  region  to  the  rear.  Only  when  there  exists  a considerable 
input  of  kinetic  energy  to  a given  mode  of  vibration  will  the  Lame 
prediction  be  exceeded.  It  is  evident  from  figure  3 that  this  con- 
dition is  not  fulfilled  in  the  breech  section  of  the  tube.  Different 
behavior  will  be  witnessed  further  along  the  bore. 
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Figure  4 - Theoretical  Displace- 
ment and  Pressure  a Short  Distance 
Beyond  the  Origin  of  Rifling. 


Figure  4 shows  the  same  set  of 
functions,  i.e.,  the  NASTRAN  and 
Lame'  predicted  displacements  and 
the  applied  ballistic  pressure,  at 
a station  located  a short  distance 
in  front  of  the  projectile  prior 
to  ignition.  The  distinctive  fea- 
ture as  compared  with  Figure  3 is 
the  time  delay  prior  to  the  arrival 
of  the  projectile.  During  this 
delay  a negative,  that  is,  an  in- 
ward radial  displacement  occurs  at 
this  location.  This  is  the  'con- 
| tractile  effect'  referred  to  in  the 
previous  paragraph,  caused  by  adja- 
cent regions  of  pressure  and  non- 
pressure. The  Lame  calculation 
cannot,  of  course,  account  for  such 
an  effect  since  all  regions  are 
pssumed  uniformly  pressurized  in 
the  Lame  problem. 


At  approximately  6.5  milliseconds  following  ignition,  the  projectile 
passes  exposing  this  portion  of  the  bore  surface  to  the  pressure  of 
the  propellant  gas.  Thus  the  time  rate  of  load  application  at  this 
station,  in  contrast  with  that  in  the  breech,  is  almost  wholly  deter- 
mined by  the  magnitude  of  the  load  velocity  at  this  time.  The  loading 
is  seen  to  be  much  more  abrupt  than  that  shown  in  Figure  3.  This  not- 
withstanding, one  notes  that  the  Lame-calculated  displacement  proves 
once  again  to  be  conservative. 


While  Figures  3 and  4 show  the  Lam^-calculated  displacement 
to  err  conservatively  at  bore  locations  nearer  the  breech,  this  is  not 
typical  of  the  behavior  at  stations  further  along  the  bore.  Figure  5 
represents  the  response  more  likely  to  be  realized  at  least  at  some 
down  bore  locations.  The  station  represented  in  this  figure  is  approx- 
imately midway  along  the  length  of  the  tube.  At  this  point  the  pro- 
jectile is  moving  with  a very  high  velocity  and  a very  abrupt  loading 
condition  results.  Provided  that  the  vibration  of  the  surface  is  in 
phase  (i.e.,  moving  outward)  with  the  pressure  application,  a relatively 
large  radial  displacement  will  result  as  is  the  case  in  Figure  5.  One 
should  note  that  only  a short  distance  forward  and  aft  of  this  location 
the  surface  movement  may  be  completely  reversed.  A sudden  application 
of  pressure  at  these  locations  would  do  negative  work  on  the  moving 
surface  and  the  oscillations  may  be  greatly  attenuated.  The  actual 
locations  along  the  tube  where  one  can  expect  overstresses  (large 
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Figure  5 - Displacement  and 


pressure  midway  along  tube 


radial  displacements  beyond  the 
Lame  prediction)  cannot  really  be 
computed  accurately  owing  to  the 
imperfect  knowledge  of  the  true 
boundary  conditions  and  material 
properties  - not  to  mention  the  in- 
ability of  any  integration  scheme 
'to  follow  the  phase  accurately  over 
a great  number  of  cycles  of  vibra- 
tion. It  is  important  nonetheless 
to  know  that  such  overstresses  are 
to  be  expected. 

As  experimental  evidence  of  the 
[success  of  the  NASTRAN  analysis,  two 
representative  test  records  from 
actual  gun  firings  are  offered  in 
*•  figures  6 and  7.  Figure  6 shows  the 
displacement  at  a location  midway 
along  the  tube  and  therefore  can  be 
compared  semi-quantitatively  with 


figure  5.  One  clearly  observes  the 

production  of  natural  oscillations  and  the  formation  of  beats  of  large 
amplitudes.  It  is  mostly  due  to  the  absence  of  damping  in  the  anal- 
ysis and  the  greater  abundance  of  modes  in  the  real  tube  into  which 


energy  may  flow,  that  the  amplitudes  in  the  test  records  are  lesser 
than  those  predicted.  It  is  notable  that  the  time-averaged  displace- 


ment estimated  from  figure  6 agrees  quite  well  with  that  of  figure  5. 
On  the  other  hand  the  excessively  large  maximum  displacement  of  figure 
5 does  not  materialize  in  figure  6.  As  previously  discussed,  this  is 


a matter  of  phase  agreement  between  the  surface  motion  and  the  applied 
pressure.  In  figure  7,  however,  an  overstress  condition  is  clearly  in 
evidence  - the  peak  Lame  value  being  indicated  for  comparison.  This 
station  is  located  nearer  the  muzzle.  One  notes  the  exceptionally 


clear  oscillation  and  beat  formations  which  continue  after  shot  ejec- 
tion which  is  estimated  to  occur  at  t = 21  milliseconds  in  the  figure. 


i 
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2.  TRANSVERSE  (BENDING)  MOTION  OF  GUN  TUBES  - CURVATURE  INDUCED  LOADS 

In  1959,  measurements  by  Gay  and  Elder  of  the  US. Army  Ballis- 
tics Laboratory  (5)  showed  that  the  muzzle  motions  of  a 90  mm  tank  gun 
at  the  time  of  shot  ejection  are  very  small  but  yet  significant  in 
explaining  a portion  of  the  error  realized  at  the  target.  Typical 
rotations  and  displacements  at  the  muzzle,  for  example,  were  of  the 
order  of  10"!  milliradians  and  10“2  inches,  respectively.  The  theory 
of  gun  tube  motion  by  which  explanations  were  sought  for  these  obser- 
vations assumed  that  the  motion  was  due  solely  to  a mass  eccentricity 
at  the  breech  which  produced  a sudden  inertial  moment  upon  recoil. 

The  theoretically  predicted  motion,  however,  was  often  much  smaller 
in  magnitude  than  that  observed.  In  the  work  that  follows  it  will  be 
shown  that  loads  resulting  from  the  initial  curvature  of  the  M-68, 

105  mm  gun  tube  can  also  produce  muzzle  motions  of  these  magnitudes 
and  should  therefore  not  be  neglected. 

Recoil  Loads 


During  the  recoil  of  a gun  tube  there  results  an  axial  load 
per  unit  length  which  is  equal  to  the  product  of  the  recoil  acceler- 
ation, a(t),  and  the  mass  per  unit  length  of  the  tube  p(x);  i.e.. 


w(x,t)  = -p(x)a(t) 


(1) 


Figure  8 - Inertial  Moment  Due  to 
Recoil . 


When  the  tube  is  curved,  this 
load  creates  a moment  at  any 
location  x,  along  the  tube. 
Referring  to  Figure  8,  the 
total  moment  is  given  by  the 
integral: 

i 

M(x,t)  = / -w(x,t) [y(x,t) 

X 


-y(x,t)]dx 


(2) 


Differentiating  twice  with  respect  to  the  space  variable  x gives  the 
resultant  transverse  load  intensity  due  to  recoil:  [ (')  = d/dx  ] 

fR(x,y*  ,y",t)  = a(t){p(x)y' (x,t)  - y"(x,t)/ P(£)d£;}  (3) 
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' Bourdon ' Load 


If  a gun  tube  is  curved,  the  bore  surface  area  b.  <nes 
asymmetrically  distributed  about  the  central  axis  owing  to  the  relative 
extension  and  contraction  of  the  material.  From  the  geometry  of  Figure 
9 one  can  verify  that  the  net  difference  in  area  is  given  by  the 
integral: 

2ir 

J f sin0+l)asin0dGdx  = -ira2y"dx  (4) 

0 

where  R"1  has  been  replaced  by  -y". 

The  applied  ballistic  pressure  there- 
fore produces  a resultant  transverse 

load  intensity:  \ 

fB(x,y",t)  = -p(x,t)7ra2y"(x,t)  (5) 

p(x,t)  is  a traveling  ballistic  pressure  ; 

function,  i.e., 

P(x,t)  = P0(t)H(5-x) 

H(z)  is  the  Heaviside  unit  step  func- 
tion and  £(t)  represents  the  distance 
traveled  by  the  projectile  along  the 
bore,  (fg  has  been  called  the  'Bourdon'  Load  because  of  its  similar- 
ity to  a straightening  Bourdon  tube  upon  pressurization.  Actually, 
the  two  effects  are  completely  different  and  the  term  'Bourdon'  load 
is  a misnomer.) 

Projectile  Loads 

If  the  projectile  is  assumed  to  be  a point  mass  mp,  travel- 
ing a bore  path  which  changes  in  time,  it  can  be  shown  (8)  that  there 
results  a transverse  load  function  containing  Coriolis,  centrifugal  and 
transverse  accelerations,  i.e.,  [ (•)  = 3/3t)  ] 

fp  = -™p[y  + 2Vy'  + V2y"  + g]6(x-C(t))  (6) 

where  <S(z)  is  the  Dirac  function  and  £(t)  = /*Vd t,  where  V is  the  pro- 
jectile velocity,  g is  the  gravitational  constant. 


s 


Figure  9 - Curved  Section  of 
Tube. 
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Equation  of  Motion 

Using  Euler-Bernoulli  beam  theory,  the  displacement  y(x,t) 
from  the  undeformed  (straight)  neutral  axis  of  the  tube  must  satisfy 
the  partial  differential  equation: 

(Ely")"  + p(x)y  -f  - (EIY") " = 0 (7) 

where  I(x)  is  the  area  cross  section  moment  of  inertia  and  E is  Young's 
Modulus  of  Elasticity.  The  last  term  on  the  right  hand  side  of  (7) 
represents  the  static  load  intensity  corresponding  to  an  initial 
deformation  Y(x). 

The  initial  conditions  are:  y(x,0)  = Y(x)  and  y(x,0)  = 0 (8) 

The  load  functions  f = fR  and  f = fg  will  be  considered  in 
turn.  It  is  to  be  noted,  however,  that  the  beam  motions  due  to  each 
load  function  may  not  be  superposed  as  the  linear  operators  involved 
vary  with  the  particular  load  function  considered. 

Taking  Y(x)  to  be  the  gravitational  deformation  of  the  tube 
prior  to  firing,  and  f = fR(x,y',y"t)  it  is  convenient  to  first  define 
y(x,t)  as  the  displacement  as  measured  from  the  initial  curve  Y(x). 

Thus  (7)  becomes: 

(Ely")"  v p(x)£  -fR(y,y  ",x,t)-f(Y',Y",x,t)  = 0 (9) 

The  initial  conditions  in  terms  of  y are  homogeneous: 

y(x,0)  = y(x,0)  = 0 (10) 

Static  and  Dynamic  Support  Conditions 

The  support  (boundary)  conditions  which  prevail  prior  to 
firing  may  not  be  appropriate  during  the  ballistic  cycle.  For  the  pur- 
pose of  calculating  the  initial  static  gravitational  deformation  Y(x), 
the  M-68  tube  is  assumed  to  be  cantilevered  from  its  two  mounting 
points  near  the  breech.  Actual  vibration  records  indicate  that  the 
mount  connections  are  far  from,  rigid,  however.  In  fact,  excellent 
agreement  between  calculated  and  experimentally  observed  natural 
frequencies  is  obtained  only  if  the  tube  is  regarded  as  virtually 
unsupported  during  the  ballistic  cycle.  The  dynamic  boundary  conditions 
are  therefore  assumed  to  be  free.  This  assumption  should  be  valid  as 
long  as  tube  displacements  are  of  the  order  of  the  support  clearances. 
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The  partial  differential  equation  (9)  was  transformed  into  a 
set  of  N - ordinary  differential  equations  in  time  via  the  Galerkin 
procedure  in  which  the  basis  functions  were  chosen  to  be  the  natural 
modes  of  vibration  of  the  unsupported  M-68  tube.  Using  once  again  the 
NASTRAN  finite  element  program,  these  mode  shapes,  W^(x)-  as  well  as 
the  initial  gravitational  deformation  function  Y(x)  - were  determined 
at  nineteen  points  along  the  tube.  Between  these  points  the  shapes 
were  interpolated  by  cubic  spline  functions.  The  approximation  for  y 
in  terms  of  the  mode  shapes  may  be  written: 

N 

y(x,t)  = l ai(t)Wi(x)  , 
i=l 


the  resulting  set  of  differential  equations  appear  as  follows: 

N 

a.  + (U.2a.  + l (A..(t)-B..(t))a.  = C.(t);  i=l  to  N (11) 

j=l  J J J 

The  initial  conditions  are  homogeneous,  i.e.,  a^(0)  = a^(0)  = 0 (12) 


ku  represents  the  natural  frequency  of  the  ith  vibration  mode.  A^ (t) 
derives  fromthe  y"  terms  in  (9)  while  B^j(t)  reflects  the  slope 
dependency,  y,  and  vanishes  for  the  case  f = fg.  C^(t)  is  due  to 
terms  in  the  initial  deformation  function  Y(x). 


The  system  represented  by  (11)  and  (12)  above,  was  solved 
numerically  using  4th  order,  variable  time  step  Runge  Kutta  integration. 
In  practice  a value  of  N = 10  gave  satisfactory  convergence.  This 
included  two  rigid  body  modes  (plane  rotation  and  translation)  and 
eight  vibration  modes  of  finite  frequency. 

RESULTS  AND  DISCUSSION 

Figures  10  and  11  represent  the  response  of  the  M-68  gun  tube 
to  recoil  and  asymmetric  ballistic  pressure  (the  'Bourdon'  effect). 

The  magnitudes  shown  in  each  of  the  figures  are  comparable  with  those 
predicted  by  Gay  and  Elder  (8)  in  connection  with  eccentric  breech 
inertia  in  the  T-139  tank  gun.  This  indicates  that  curvature-induced 
loads  should  be  included  in  any  theory  of  gun  tube  motion  during  firing. 

A fundamental  difference  is  apparent  between  figures  10  and 
11,  and  is  especially  obvious  in  comparing  the  time  histories  of  the 
tube  shapes,  i.e.,  figures  10a  and  11a.  Whereas  in  10a  the  motion  is 
well  developed  along  the  entire  tube  length,  the  motion  in  11a  is  much 
more  wavelike  and  muzzle  displacements  remain  comparatively  small 
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through  shot  ejection.  This  difference  is  due  to  the  fact  that  the 
recoil  inertia  load  acts  instantaneous,  y*  over  the  entire  tube  length 
whereas  development  of  the  'Bourdon'  load  proceeds  at  projectile  velo- 
city. Moreover,  the  recoil  loading  consists  of  two  parts  (cf.  equa- 
tion 3).  If  evaluated  near  time  zero,  these  two  parts  consist  of  a 
downward  load  proportional  to  the  curvature  (y")  and  therefore  strong- 
er near  the  supports,  and  an  upward  load  proportional  to  the  slope  (y') 
more  intense  near  the  muzzle.  The  two  act  in  unison  to  encourage  a 
rotation  of  the  entire  tube.  The  large  mass  at  the  breech,  however, 
tends  to  anchor  that  end  of  the  tube.  The  result  is  the  deformation 
shown  in  figure  10a.  The  Bourdon  load,  on  the  other  hand,  consists  of 
only  one  part  (cf.  equation  5).  This  part,  being  proportional  to 
curvature  is  greatest  near  the  supports  and  trails  off  to  practically 
zero  value  by  mid-length  of  the  tube.  Muzzle  displacement  must  there- 
fore await  the  arrival  of  a disturbance  originating  from  the  support 
end  of  the  tube. 

A similarity  is  also  apparent  comparing  figures  10  and  11. 

The  time  histories  of  the  slopes  at  the  muzzle  - probably  the  most 
important  of  all  the  response  curves  - show  a definite  similarity  in 
that  little  muzzle  rotation  occurs  during  the  first  half  of  the  bal- 
listic cycle.  This  delay  represents  the  propagation  time  for  a dis- 
turbance established  near  the  supports  (where  y"  is  largest)  to  reach 
the  muzzle. 

CONCLUSION 

To  completely  account  for  the  motions  of  a gun  tube  prior  to 
shot  ejection  several  effects  must  be  considered,  among  them  the  effects 
of  tube  curvature.  Of  particular  significance  is  the  muzzle  slope  at 
the  time  the  projectile  arrives.  Slope  values  are  particularly  sen- 
sitive to  sudden  disturbances  which  originate  near  the  support  or 
breech  end  of  the  tube  and  propagate  in  wavelike  fashion  reaching  the 
muzzle  prior  to  shot  ejection.  This  is  exemplified  in  the  case  of 
the  Bourdon  load  and  strongly  suggests  an  opportunity  for  further 
study  as  the  supports  themselves  may  create  strong  forces  as  clearances 
are  abruptly  taken  up  during  the  initial  phases  of  recoil.  The  result 
may  significantly  strengthen  rotations  of  the  muzzle  at  shot  ejection. 


•The  tube  is  assumed  axially  rigid  in  the  analysis. 
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INCHES  FROM  REAR  FACE  OF  TUBE  TIME  (MILLISECONDS) 

Figure  11a  - Displacement  of  Tube  Figure  lib  - Muzzle  Displacement 
Axis  from  Initial  Position. 


TIME (MILLISECONDS)  TIME (MILLI5EC0N0S) 


Figure  11c  - Muzzle  Slope  Figure  lid  - Transverse  Velocity 

(Rotation).  of  Muzzle. 


Figure  11  - Transverse  Response  of  the  M-68  Gun  Tube  to  Traveling 
Bourdon  Pressure. 
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